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ABSTRACT 

Long-duration gamma-ray bursts (GRBs) provide a premier tool for studying high-redshift star-forming 
galaxies thanks to their extreme brightness and association with massive stars. Here we use GRBs to study 
the galaxy stellar mass-metallicity (M*-Z) relation at z ~ 3 - 5, where conventional direct metallicity measure- 
ments are extremely challenging. We use the interstellar medium metallicities of long GRB hosts derived from 
afterglow absorption spectroscopy (Z ^ 0.01 - 1 Z©), in conjunction with host galaxy stellar masses determined 
from deep Spitzer 3.6 jim observations of 20 GRB hosts. We detect about 1/4 of the hosts with MabW ~ -21 .5 
to -22.5 mag, and place a limit of MabW ^ -19 mag on the remaining hosts from a stacking analysis. Using 
these observations, we present the first rest-frame optical luminosity distribution of long GRB hosts at z ^ 3 
and find that it is similar to the distribution of long GRB hosts at z ~ 1 . In comparison to Lyman-break galaxies 
at the same redshift, GRB hosts are generally fainter, but the sample is too small to rule out an overall sim- 
ilar luminosity function. On the other hand, the GRB hosts appear to be more luminous than the population 
of Lyman-alpha emitters at z ~ 3-4. Using a conservative range of mass-to-light ratios for simple stellar 
populations (with ages of 70 Myr to ~ 2 Gyr), we infer the host stellar masses and present mass-metallicity 
measurements at z ~ 3 -5 ((z) ~ 3.5). We find that the detected GRB hosts, with ^ 2 x 10^^ M©, display 
a wide range of metallicities, but that the mean metallicity at this mass scale, Z ^ 0.1 Z©, is lower than mea- 
surements at z ^ 3. Combined with stacking of the non-detected hosts with < 3 x 10^ and Z < 0.03 
Z0, we find evidence for the existence of an M^-Z relation at z ^ 3.5 and continued evolution of this relation 
to systematically lower metallicities from z ~ 2. 
Subject headings: 



1. INTRODUCTION 

The simple ' ' closed -box" model of galaxy evolution 
(iTalbo t & AmettI Il97l|) predicts a correlation between the 
stellar mass and the gas-phase metallicity of a galaxy (the 
M^-Z relation), under the assumptions of no gas inflows 
or outflows, a constant yield of metals, an invariant stel- 
lar initial mass function (IME), and instantaneous mixing of 
newly- synthesized metals back into the interstellar medium 
(ISM). In reality, this simple picture is complicated by 
the fact that galaxies accrete low-metallicity gas from the 
intergalactic medium (IGM) and lose metal-enriched gas 
through galactic- scale winds or by depletion on to dust. 
In addition to these processes, the M^-Z relation may also 
be modified by a mass-dependent star-formation efficiency 
jjuneau et al.l i 2005: Ee ulner et al.l 120051: lErancesc hini et alj 
l2006l: lAsari et al.. ■2007h. and poss ibly an environmental- 
dependent IME (iKoppen etalJl2007h . Thus, the M*-Z rela- 
tion and its evolution with redshift provide insight into the 
physical processes that shape galaxy formation and evolution 
across cosmic time. 

Given the importance of this relation it has been the focus 
of several exte nsive studies out to z ~ 3. In the local uni- 
verse (z ^ O.D. lTremontiet al .' ("2004) studied 53,400 galaxies 
from the Sloan Digital Sky Survey and found a tight corre- 
lation (±0.1 dex) between stellar mass and metallicity over 
a range of ^ ^q8.5 -ii.5 order o f magnitude 

in metallicity (see also iKewley & Ellisonll2008|) . They con- 
cluded that the observed correlation is best explained by the 
influence of metal-enriched outflows, with a larger metal loss 

^ Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cam- 
bridge, MA 02138 

^ Spitzer Science Center, California Institute of Technology, Pasadena, CA 
91125 



in lower mass galaxies. Stu dies at z 7^ 1 - 3 (ISavaglio et al 



2005; Erbetal. 2006; Maioli no et all 120081: iMannucci et al.. 

■2009: Jabran Zahid et al 20 id) found that the M*-Z relation 

evolves by about 0.8 dex fro m z ^ 3 to the presen t, while 
keeping the same overall trend. 'Savagl io et al.l (12005) argued 
that the redshift evolution to z ~ 0.7 can be reproduced in the 
simple closed-box model with the assumption that the typical 
timescale for star f ormation is longer in lower mass galaxies. 
On the other hand, iJabran Zahid et al.l (l2010l) found that the 
M^-Z relation evolves only below ^ 10^^-^ to z ~ 0.8, and 
argued that unlike in the local universe the effective yield de- 
creases with larger mass and that a closed-box model cannot 
explain the evolution. They further argued that outflows play 
a minor role, and proposed that a rising star forma tion effi- 
ciency with large mass may be the dominant effect. Erb et al.l 
(l2006h argued that at z ~ 2.3 the primary driving mechanism 
for the -Z relation and its evolution is the expected increase 
in metallicity as star formation leads to a reduced gas fraction, 
and that outflows affect galaxies at all mass scales. 

Beyond z ~ 2.3 there are o nly a few m easurements of 
galaxy metallicities and masses. IMannucci et al. ( 2009) stud- 
ied 10 Lyman break galaxies (LBGs) at z ~ 3.1 and found 
continued downward evolution of the M^-Z relation, and a de- 
creasing effective yield with larger stellar mass. They argued 
that gas infall plays the domina nt role in the M^-Z rel ation, 
and that outflows are not needed. iMaiolino et al.l (l2008h stud- 
ied 9 LBGs at z ~ 3.5 and found a decline in the mean metal- 
licity at a stellar mass scale of 1.4 x 10^^ compared 
to z ~ 2.3, with a possible steepening of the M^-Z relation 
relative to lower redshifts. 

Tracing the M^-Z relation and its evolution to even earlier 
times will provide insight into the earliest epochs of galaxy 
evolution, while allowing us to probe the relative importance 
of the various galactic- scale phenomena proposed at z ^ 3. 



2 



Laskar et al. 



Although initial studies of LBGs at z ^ 3.5 are now available 
(iMaiolino et al .120081: iMannucci et al.ll2QQ9l) . these studies are 
challenging because the nebular emission lines required for 
robust metallicity measurements^ (e.g., Ha, H/3, N IIA6583, 
O IIIAA4959, 5007, O IIAA3726, 3729) shift into the near- and 
mid-IR, where existing spectrographs have reduced sensitiv- 
ity compared to the optical band. This is further complicated 
by the rapid dimming of galaxies at higher redshift such that 
only the most luminous LBGs are amenable to spectroscopy. 

An alternative way to determine metallicities_atz> 3 (and 
in principle at z ^ 10 and beyond; Salvaterra et al.l I2 QQ9 I; 
iTanvir et al.l l2009|) is absorption spectroscopy of gamma-ray 
burst (ORB) optical/near-IR afterglows. Long-duration GRBs 
are kn own to be associated with the deaths of massive stars 
(e.g., Woosley & B loomll2006 ). and therefore with sites of ac- 
tive star formation. The large optical luminosities of ORB 
afterglow s (easily exceeding 20 mag for several hours even 
at z ~ 8; iTanvir et aP |2009|) . and their intrinsic featureless 
spectra, provide a unique way to measure interstellar medium 
(ISM) metallicities for galaxies at z ^ 2 from rest- frame ultra- 
violet metal absorption lines and Lya absorption. Since the 
afterglows are significantly brighter than the underlying host 
galaxies, this technique allows us to measure metallicities in- 
dependent of the galaxy brightness. Moreover, since long 
GRB progenitors reside in star forming environment within 
their hosts, their sight-lines probe the warm ISM and H II re- 
gions that give rise to the (rest- frame optical) nebular emis- 
sion lines that are used for metallicity measurements at z < 3. 
This approach has now been exploited at least to z ^ 5 using 
optical spectra (e.g.. Berger et al.l2006HPro chaska e t al.l2007l: 
[pynbo et al. 2009), and with near-IR spectrographs it can be 
implemented to z ~ 20. 

Naturally, to explore the M^-Z relation at z ^ 3 we also re- 
quire a determination of the GRB host galaxy stellar masses, 
and hence follow-up infrared observations with the Spitzer 
Space Telescope to probe the rest-frame optical luminosity. 
Here, we present the first large set of Spitzer observations 
for GRB host galaxies at z ~ 3 - 5, and combine the inferred 
masses with measured metallicities to explore the M^-Z re- 
lation beyond z ~ 3. Since deep Spitzer/IRAC images are 
generally confusion-limited, our use of GRB afterglows pro- 
vides an additional boon — they accurately pinpoint the lo- 
cation of the host galaxies (to ~ O'M), thereby allowing for 
accurate galaxy identifications^. The plan of the paper is as 
follows. We present the Spitzer observations, analysis, pho- 
tometry, and metallicity data in ^ In ^3.11 we present the 
first rest-frame optical luminosity distribution of GRB hosts 
at z ^ 3 and compare it to both z ~ 1 GRB hosts and field 
galaxy samples at z ~ 3. We derive the mass distribution in 
^3.21 Finally, in ^we combine the mass and metallicity mea- 
surements to place the first points on the M^-Z diagram at 
z ^ 3. We explore the implications of our results and future 
prospects in ^ 

^ For example, 'Mannucci et al." ('2009') use the R23 diagnostic to determine 
metallicities for their z ~ 3 galaxy sample, but this relation is known to be 
double- valued. They attempt to discriminate between the low- and high- 
metaUicity branches using the O IIIA5007/O IIA3727 ratio. 

This can be contrasted with the potential use of quasar intervening ab- 
sorption systems for studies of the M*-Z relation, since the galaxy coun- 
terparts of the absorbers are offset on the sky by ~ few arcseconds from 
the quasar position. Spectroscopic confirmation is therefore required to de- 
termine the correct counterpart, negating the advantage of absorption spec- 
troscopy. Furthermore, even if the galaxy counterparts could be identi- 
fied, observations with Spitzer's large point- spread-function against the much 
brighter quasar glare are essentially impossible. 



2. GRB SAMPLE AND DATA ANALYSIS 

We obtained deep observations of all 35 long GRB host 
galaxies with spectroscopic redshifts in the range z ~ 2-5.8 
available as of November 2006 u sing the 3.6 /im b and of the 
Infra-Red Array Camera (IRAC; iFazio et al1l2004l) on-board 
the Spitzer Space Telescope. Here we investigate the proper- 
ties of GRB hosts in the redshift range z ^ 3 - 5.6 (observed 
for about 2 hours each); targets at redshifts z ^ 2-3 are treated 
elsewhere (Chary et al. 2011, in prep). For the objects in 
this paper, the effective wavelength of the IRAC 3.6 /im band 
probes the rest-frame spectral energy distribution (SED) red- 
ward of about 5500 A (Figure [TJ and therefore provides a ro- 
bust measure of the stellar mass. 

We processed the Sp itzer data using the standard mopex 
(iMakovoz et al.l l2006h software package to generate mo- 
saics for each target. The mopex package detects and 
removes cosmic rays and moving objects before drizzling 
(iFruchter & Hookl 12002). performing background equaliza- 
tion, and applying distortion corrections. For our coverage 
and dither pattern we find that an output pixel scale of 0'^4 
and a drizzling parameter of 0.7 provide the best combination 
of improvement in the point- spread-function (PSF) with min- 
imal degradation of the signal-to-noise ratio. We set all other 
parameters in MOPEX to their recommended defaults. 

2.1. Astrometry 

We used optical afterglow images to perform relative as- 
trometry on the Spitzer mosaics and to locate the GRB hosts. 
The median root-mean- square residual of the astrometric ties 
is about 0'M2, corresponding to about one-tenth of the Spitzer 
PSF at 3.6 /im. This is the dominant source of uncer- 
tainty in the astrometry, since the optical afterglow detec- 
tions themselves are mostly of high signal-to-noise. In only 
the two cases (GRBs 050502 and 050814), where no af- 
terglow images were available, we performed absolute as- 
trometry based on SDSS and 2MASS using the afterglow 
coordinates as reported in t he GCN circulars (Jensen et al.l 
2005; Blak e &Bloomll2005h . We detect one of these hosts 
(GRB 050814) in our Spitzer follow-up. The Spitzer images 
for the detected hosts are presented in Figure |2] while non- 
detections are presented in Figure |3] 

2.2. Photometry 

At the depth of our observations, Spitzer images are 
confusion-limited for faint sources. As a result, in several 
cases the region around the expected location of the GRB host 
is contaminated by light from nearby stars or galaxies. Prior to 
perf orming photomet ry, we used the GALFIT software pack- 
age (iPeng et al.ll2002 h to model and subtract these neighbor- 
ing sources. For this purpose, we used multiple point sources 
to determin e the mosaic P SF with the IDL STARFINDER rou- 
tines ( Diolaiti et al.ll2000|) . The accuracy of the generated PSF 
was evaluated by fitting and subtracting point sources at vari- 
ous locations on the mosaic. To remove neighboring sources 
around the expected location of the hosts, we used GALFIT 
with point-source, Gaussian, or Sersic models as appropriate 
in order to achieve the lowest level of residuals. In two cases 
(GRBs 061222B and 050505), the expected location of the 
host based on the afterglow astrometry fell on the diffraction 
spike of a saturated star. Since we cannot model the PSF at the 
required level of accuracy to robustly subtract these diffrac- 
tion spikes, we do not consider these two sources in the sub- 
sequent analysis. 
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For the remaining 18 targets, we searched within 1" .5 of 
the afterglow centroid (corresponding to ^ 10 kpc at z~ 3.5) 
and detected 5 GRB hosts at 3.6 /am (GRBs 050319, 050814, 
060707, 060210, and 060926). None of these hosts were de- 
tected in the simultaneously-observed 5.8 fim IRAC band, 
which has substantially worse sensitivity. In the two cases 
(GRBs 060926 and 060210) where a nearby source was sub- 
tracted prior to photometry, we found (based on the level of 
residuals) that the flux we associate with the GRB host cannot 
be explained by modelling it as part of the subtracted source. 
For GRB 050908, a visual inspection reveals a coincident flux 
excess, but photometry indicates that it is consistent with a 
noise fluctuation. For GRB 0600206, the source V\5 to 
the West of the afterglow centroid is an unrelated foreground 
object. 

To estimate the probability that one or more of the detected 
sources are chance superpositions, we ran STARFlNDER's 
source-detection routines on 3^4 x 3^4 pixel cutouts of the 
field around our targets and searched for sources down to 5 a 
(the significance level of the faintest detection) using the PSF 
generated from the corresponding images. Based on the mean 
number o f sources dete cted at different thresholds and follow- 
ing Bloo m et a we assign a false detection probabil- 
ity given by P^c = 1 - ^^"^ to each of our detections (Ta- 
ble [T]). Here Pec is the probability of chance coincidence, R 
is the aperture radius and T>(u) is the number of sources per 
unit area down down to the threshold, u = ka. The probability 
that all of our detections are chance coincidences is negligi- 
ble (10~^), while the probability that none of the targets are 
chance superpositions is 72%. Thus, whereas it is possible 
for one or two of our detections to be chance superpositions, 
it is highly unlikely to be the case for all. 

We use the funtools package to perform aperture pho- 
tometry on our detections by placing apertures of 2 native 
IRAC pixels (2.45^0 in radius and background annuli of 2 - 6 
native pixels (2.45-7.34^0 in radius centered on the detected 
sources. We choose these values since they allow us to ap- 
ply standard IRAC aperture corrections^, which are relevant 
for the expected compact sizes of galaxies at z > 3. In the 
cases where this choice of radii caused nearby objects above 
the 3(7 level to fall within either the aperture or the back- 
ground annulus, we mask them out and correct for the lost 
flux (in both the aperture and annulus) by determining our 
own aperture correction using mosaics of the IRAC calibra- 
tion star HD 18 12095, prepared with identical parameters as 
for our targets. 

We determine uncertainties on our measured flux densi- 
ties using the uncertainty mosaics created by MOPEX. We 
carry out aperture photometry in an identical fashion on the 
(squared) uncertainty images as for the source images them- 
selves, including aperture corrections as described above. In 
addition, we account for correlated noise due to the drizzling 
process by incorporating an estimate for it in the flux density 
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uncertamty": 



A72 



B 



(1) 



where oSyc is the variance of the source flux density, A is the 
aperture correction, Fcon is the effective number of pixels over 

^ See section 4.10 of the IRAC Instrument Handbook. 
^ See http://wise2.ipac . caltech. edu/staff/fmasci/ ApPhotUn cert_corr.pdf | 
for a derivation. 



which noise is correlated in the mosaic, A^a and A^b are the 
number of pixels in the aperture and background region, re- 
spectively, and CT/^A and ct/^ are the uncertainty of the flux 
density in the /* pixel of the aperture and background regions, 
respectively. Since our final mosaic has 0"A pixels, whereas 
the native detector pixels are 1''.22 on a side, noise will be 
correlated over about 3 pixels in our images (the exact corre- 
lation function will depend on the drizzling parameter as well, 
for which we used a value of 0.7). As a conservative estimate, 
we take Fcorr to be 3. 

The resulting flux densities and upper limits are listed in 
Table [U The detections range from about 0.55 to 1.65 /iJy, 
while the typical upper limit is about 0.25 /iJy (3a). We 
list the Spitzer 3.6 jim result s for five additiona l GRB hosts 
at z ^, 3 from the literature (Berg er et al.ll2007i IChary et aP 
[2007l:lchen et"ID[2010) in Table[2l 

2.3. Stacking 

To assess the typical flux density of the non-detected hosts 
we carry out a stacking analysis with 11 of the 13 non- 
detections that have accurate relative astrometry. We exclude 
GRB 050502, for which we only have absolute astrometry, 
and GRB 060927, for which the relative astrometry is poor 
due to a low signal-to-noise detection of the afterglow. The re- 
maining 1 1 targets are located at a median redshift of z = 3.4. 
We first perform sub-pixel shifts on the Spitzer mosaics to 
bring the expected location of each host (based on the af- 
terglow centroid) to the center of a mosaic pixel using the 
IRAF task imshift. We then average 51x51 pixel sections 
from each image centered on the expected location of the host, 
weighted by the inverse of the corresponding variance maps, 
after masking out the bright (> 5cr sources). The resulting 
stacked image (Figure |4]) does not show a detection, and we 
place a limit on the mean flux density of these 11 hosts of 
< 80 nJy (3cr). 

2.4. Metallicities 

Absorption spectra of GRB afterglows present a unique op- 
portunity to measure ISM abundances of galaxies at z ^ 3, 
where current spectroscopic sensitivity limits are inadequate 
for measuring metal abundances. A typical optical afterglow 
spectrum exhibits a wide range of ISM absorption features 
due to rest-frame UV transitions of low- and high-ionization 
metal species, which allow a direct determination of the col- 
umn density of these elements along the GRB line of sight 
through the host galaxy. Combined with a determination of 
the neutral hydrogen column density via the L yg line, it is 
possible to determine the ISM abundances (e.g . , iBerger et"an 
l2006l:IProchaska et al.ll2Q07l: iFvnbo et al."2009'). 

Of the ions typically present in an afterglow spectrum, 
many are due to refractory elements and therefore depletion 
onto dust precludes their use as robust abundance indicators 
(they can be used to place a lower limit on the metallicity). 
In this work we use S II, when available, as a measure of the 
metallicity, primarily since sulfur is not strongly depleted onto 
dust. Furthermore, the S IIAA1250.6, 1253.8 transitions have 
low oscillator strengths, and the corresponding lines are more 
likely to be unsaturated. 

In Table [3] we present a compilation of metallicities for 
our GRB host sample, including the spectral line used, the 
neutral hydrogen column density, and the column density 
of the metal ion computed using the measured equivalent 
width of the transition. All except one of these values are 
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taken from the literature (i Hiorth et"a Dl2003t 'Vreeswiik et alJ 
2QQ4I: ISav aglio 2006; Ferre ro et al.l 120 09: Fynbo et al. 2009; 
Ledoux et al.,.2009.: .Thone et alJ 20 IQ) and have bee n placed 
on the solar abundance scale of lAsplund et alJ (l2005h . Where 
an S II line was not detected, we use the errors reported on the 
equivalent widths of transitions in the vicinity of rest-frame 
1250 A to place 3cr upper limits on the metallicity. We also 
place lower limits on the metal licity u s ing Si II, Si IV and 
C II detections reported by Fynb o et aH (l2009b . Finally, for 
five GRBs in our sample, metal line equivalent widths are not 
available and we do not consider them in our mass-metallicity 
analysis (® . The same quantities for the five additional GRB 
hosts at z ^ 3 from the literature are listed in Table |4] 

3. OPTICAL LUMINOSITIES AND STELLAR MASSES OF LONG 
GRB HOSTS AT Z 3 -5 

Before we address the M*-Z relation itself, we explore the 
rest- frame optical properties of our z^3-5 long GRB sample 
since this analysis has not been performed previously. We 
compare our sample to long GRB hosts at z ^ 2 to explore any 
evolution in host properties, as well as to field galaxy samples 
at similar redshifts to place the long-duration GRB hosts in a 
broader context. 

3.1. Luminosity Distribution 

Since the Spitzer 3.6 jim band probes different parts of the 
host SEDs at different redshifts (Figure [T}, we must correct 
the inferred luminosities to a common rest-frame wavelength 
for a meaningful comparison (K-correction). Doing so re- 
quires knowledge of the host SED, which we do not have for 
our targets. We therefore employ evolutionary single stellar 
population (SSP) models with a single burst of star forma- 
tion (e-folding tinie, r = 0) to determine the K-corrections. 
iLeibler & Bergerl (|2010|) recently performed stellar popula- 
tion modeling of 23 long GRB hosts at z ^ 0. 3 - 1. 6 us- 
ing multi-band photometry from 'Savaglio et alj (l2009h and 
the evolutionary models of Maraston (2005). They deter- 
mined a median stellar p opulation age of 10~^-^^^-^ Gyr (see 
also S avaglio et al.l 12009 ). Taking this age range into ac- 
count, along with a Salpeter IMF, and a metallicity range of 
0.05-0.5 Z©, we find that the flux density of the SSP mod- 
els in the rest-frame 0.4-2 (im range is well- approximated 
by a power law. Fx oc X^, with /3 = -2.3 ± 0.2. This assumes 
negligible extinction, which is appropriate for most GRB host 
galaxies ( Savaglio et al.ll2009|) . Using this result as an esti- 
mate of the intrinsic spectrum of GRB hosts at z > 3, we de- 
termine the host absolute magnitudes including K-correction 
as follows: 



MAB(Ao) = mAB-51og 



-2.5(/3 + 2)log 
+2.51og(l+z), 



10 pc 

(l+z)Ao 



A 



(2) 



where Mab is the absolute magnitude at a rest wavelength, Aq, 
to which the K-correction is performed, di is the luminosity 
distance, and A is the effective wavelength of the observed 
band (3.55 fim for our Spitzer observations). To minimize the 
K-corrections for our sample we select a nominal rest-frame 
wavelength, Aq = 7800 A (Figure [B, roughly corresponding 
to the /-band. The resulting mean K-correction on our pho- 
tometry is about -1.7 mag (dominated by the last term on the 



right-hand- side of Equation O, with a standard deviation of 
about 0.2 mag. 

The resulting luminosity distribution for our sample is 
shown in Figure [51 and the inferred absolute magnitudes 
are listed in Table \T\ Also included are the 5 host galax- 
ies from the literature (their values are listed in Table |2}. 
In the comparison to low redshift GRB hosts and to field 
galaxies provided below we treat separately our uniform sam- 
ple, and the combined sample that includes the 5 hosts from 
the literature. For our detected hosts we find a range of 
MAfi(780nm) ^ -21.5 to -22.5 mag, while the limits are 
typically MAfi(780nm) > -20.4 mag. The stack limit corre- 
sponds to MAfi(780nm) > -19.1 mag. With less than a 50% 
detection fraction it is not possible to robustly estimate the 
median luminosity of our sample, but the formal 3 a upper 
limit is about -21.5 mag. The stack limit, however, suggests 
that a more reasonable upper bound on the median is about 
MAfi(780nm) > -19.1 mag. The addition of the 5 hosts from 
the literature does not change this result, wit h the exception o f 
the highly luminous host of GRB 080607 ( Chen et al.ll20TQh . 
which was specifically targeted due to evidence of large after- 
glow extinction. 

3.1.1. Comparison with Long GRB Hosts at 

To compare the resulting optical luminosity distribution to 
GRB hosts at low redshift, we obtain a comparison sam- 
ple from ground-based JHK and Hubble Space Teles cope 
814 nm photometry re ported in lCastro Ceron et aD(l2008l) and 
ISavaglio et aP (l2009l) . and from our own GRB host follow- 
up studies. The comparison sample has a redshift range of 
z ~ 0.01 - 2 with a median of z ~ 0.86. For each host we 
select the band that corresponds most closely to a rest frame 
wavelength of 7800 A to minimize the K-corrections relative 
to our Spitzer sample. The mean and standard deviation of the 
K-corrections are -0.6 and 0.3 mag, respectively. 

The luminosity distribution for the low redshift sample is 
shown in Figure O There is clear overlap between the two 
GRB host samples at the bright end, but due to the lower red- 
shifts of the comparison sample, its luminosity distribution 
extends to much fainter levels (reaching ^ -16.5 mag), with a 
median of about -20.1;^q 4 mag (95% confidence range). This 
is consistent with the upper bound on the median luminos- 
ity of our z ~ 3 - 5 sample. Indeed, a log-rank test including 
the individual non-detections indicates that the hypothesis that 
two samples are drawn from the same underlying populations 
has a p- value of 0.65. Similarly, the fraction of detected hosts 
above our threshold of Mas (7 80 nm) ^ -21.5 mag is about 
1/4 for both samples. Finally, our stack non-detection level 
of Mas (7 80 nm) > -19.1 mag is consistent with the median 
of the z ~ 1 GRB sample. Thus, we find no evidence for sig- 
nificant evolution in the optical luminosity function of GRB 
hosts from z 1 to ~ 4. We note that dividing the compari- 
son sample into z < 1 and z > 1 subsets does not change this 
result. 

3.1.2. Comparison with Lyman-Break Galaxies and Lyman-alpha 

Emitters 

To assess whether the luminosities we find for GRB hosts 
at z ^ 3 are typical of field galaxies, we next compare the re- 
sulting luminosity distribution with other galaxy samples at 
similar redshifts: Ly man-break galaxies (LBGs) and Lyman- 
alpha emitters (LAEs). For the LBG comparison we use 
Spi tzer 3.6 iim photometry in the GOODS - N field reported 
by iReddv et al.1 (l2006l) : iMagdis et aH (l2010h . based on deep 
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(~ 95-hour) Spitzer observations. We also inc lude a sample 
of 72 s tar-forming galaxies atz~2.3±0.3 from lShapley et alJ 
(2005') selected based on their rest- frame UV brightness. The 
resulting luminosity distributions for these samples are shown 
in Figure [6l and their summary statistics are listed in Table [3 
K-corrections for the SED shape (third term on the RHS of 
equation O have not been applied to any of the samples, al- 
though the difference between the K-corrections should be 
minor (^0.2 mag) and would not modify the shape of the 
distributions. 

Our GRB host sample is clearly missing the luminous tail 
of LBGs atMAB ^ -23 mag, which accounts for ^ 20% of the 
comparison samples. Even if we include the luminous host of 
GRB 080607, it accounts for only 4% of the GRB host sam- 
ple observed with Spitzer. The median absolute magnitude 
of the LBG sample is -22.0!^;^ mag for R eddy^t al. (2006) 
and -21. 9!^:^ mag for lMagdis et'aP (l2010h (95% confidence 
ranges), brighter than the 3 a upper limit for the GRB sample 
(-21.5 mag, even without K-corrections for the SED shape) 
and our stack limit. A log-rank test yields a p- value of 0.23 
and 0.38 that the GRB host sample is drawn from the same 
population as the p aren t population of the LBG samples of 
iReddv et al.] (l2006l) and'Ma gdis etal.1 (l201Qh . respectively. It 
is important to note, however, that since the LBG sample is 
flux limited (based on the initial optical selection and spec- 
troscopic confirmation), whereas the GRB host sample is not, 
we cannot simply use the fractional detections of LBGs as an 
indication of the overall luminosity function. 

To further assess whether the dearth of GRB hosts with 
^AB ^ -23 mag in our sample is significant, we instead need 
to integrate the rest-frame optical luminosity function. This 
will allow us to assess the expected fraction of GRB hosts 
with Mar 5-. -23 mag compare d to our threshold of about -21 
mag. iMarchesini et al.l (l2007h calculated a Schechter fit to 
the V-band luminosity function at z ~ 3 and found M^g(y) = 
-22.77 ± 0.22 mag and faint-end slope, a = -1 . 12 ± 0.24. Ap- 
plying a K-correction from V to / band using Fx ex A~^^^^^ 
(gTB we find M^b(/) = -22.65 ± 0.30 mag. Assuming that a 
is the same in the Lband as in the V-band, we find that about 
10% of our sample (or about 2 hosts) should have Mab ^ -23 
mag if the GRB hosts are drawn from the field galaxy popu- 
lation^. Since this small number is fully consistent with zero 
detections, we cannot rule out the hypothesis that GRB hosts 
are drawn from the general LBG population. 

For the comparison to Lyman-alpha emitters we use the 
sample of lOno et aP |2010), which includes 205 LAEs at 
z ~ 3.1 and 67 LAEs at z ~ 3.7 with multi-band photometry. 
These authors find 1 1 detections at 3.6 fim (5 at z ~ 3.1 and 6 
at z ~ 3.7), corresponding to a detected fraction of only 4%. 
From a stacking analysis of the non-detected LAEs they deter- 
mine (Mab) = -20.8 mag at z ~ 3.1 and (Mab) = -21.1 mag 
at z ~ 3.7. A comparison between the GRB hosts and LAEs 
at z ~ 3.1 - 3.7 is shown in Figure |7l The luminosity dis- 
tributions of the detected LAEs and GRB hosts appear to be 
consistent. However, including the LAE non-detections and 
carrying out a log-rank test we find a negligible probability 
that the GRB host sample is drawn from the same population 
as the LAE sample at z 3.4 since our sample has a much 

^ We verify this approach by calculating the expected fraction of galaxies 
with Mab % -23 mag for our various comparison samples. We find that rela- 
tive to their typical threshold absolute magnitude of about -21.5 to -22 mag 
(Figure |6), this expected fraction is about 20%, which is in good agreement 
with the observed fraction. 



higher detection fraction than the LAEs. This is an interest- 
ing result since seve ral GRB hosts have been previously de- 
tected as LAEs (e.g. jFvnbo et al.l2QQ2Ll2003l:lJakobsson et al.l 
l2005h . It suggests that GRBs select the more luminous end of 
the LAE luminosity distribution. 

To conclude, the comparisons to LBGs and LAEs suggest 
that GRB hosts at z ^ 3 are currently missing the bright end 
of the LBG luminosity distribution (with the exception of the 
host of GRB 080607), but that this may be due to the small 
sample size. On the other hand, GRB hosts sample the high 
end of the LAE luminosity distribution. 

3.2. Stellar Mass Distribution 

We next turn to a derivation of the stellar masses of our 
GRB host sample. Computing stellar masses from observed 
luminosities in a given wave-band requires knowledge of the 
mass-to-light ratio and hence the stellar population age and 
metallicity. When multi-band photometry is available, model- 
ing of the spectral energy distribution (SED) using stellar pop- 
ulation synthesis models can used to determine stellar masses, 
provided that a single stellar population is assumed. When 
multi-band photometry is not available, the resulting uncer- 
tainty in the mass-to-lig ht ratio (e.g., at ^ 1 jim) is about an 
order of magnitude (e.g. jMagdis et al.lbOlOh . 

Here, since we lack broad-band photometry, we determine 
a range of mass-to-light ratios for each galaxy in the observed 
3.6 jim band using a wide range of population ages and the 
single stellar population models of iMarastonl (l2QQ5h with a 
Salpeter IMF. We assume an instantaneous burst of star for- 
mation (r = 0). As expected, the 3.6 jim mass-to-light ratio 
for these models increases with stellar population age beyond 
'^10 Myr. This is shown in Figure [8] where the ratio of the 
stellar mass to the observed flux density at 3.6 fim is plot- 
ted as a function of age and redshift. The more traditional 
mass-to-light ratio in solar units (in the rest-frame /-band) is 
also plotted for comparison. The upper bound on the mass-to- 
light ratio is achieved by setting the stellar population age to 
the age of the universe at each host redshift (^ L8 Gyr at the 
median redshift of our sample). We stress that this leads to a 
very conservative maximum mass for each host galaxy since 
studies of LBGs and LAEs indicate typical p opula tion ages 
of ^ 0.1-0 6 Gyr (Shaplev et al. 2005|; iReddv et al. 1120061: 

iMagdis et al.|[2010. : Qno et al 201 0). For a more typical age 

we adop t the median age for long GRB hosts at z ~ 1 of about 
70 Myr (iLeibler & Berger "2010V The variation in mass-to- 
light ratio between these age values is about an order of mag- 
nitude, as expected from other galaxy studies. The 70 Myr 
and maximum mass-to-light ratios are listed in Table [T] for 
our sample, and in Table [2] for the 5 hosts from the literature. 

We test the effect of metallicity on the mass-to-light ratio 
by considering population synthesis models at 0.02 and 0.5 
Z0, which cover the typical range of GRB host galaxy metal- 
licities. The resulting variation in mass-to-light ratio is only 
~ 15%, with the lower metallicity models typically yielding 
systematically smaller values (although this effect is redshift- 
dependent). This is a much smaller effect than the uncertainty 
due to the unknown stellar population age. In the following, 
we adopt the mass-to-light ratios for a metallicity of 0.02 Z©. 

The inferred masses of our GRB host sample are plot- 
ted as a function of redshift in Figure [9] and are listed in 
Table [T] The maximal masses inferred for our sample are 
(2.5 -5.8) X 10^^ M0, while the typical (maximal) upper lim- 
its are < 9 X 10^ M©. The masses inferred for a 70 Myr old 
population are about (0.6-1.4) x 10^^ M©, with typical up- 
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per limits of < 2 x 10^ Mq. The mass limit from the stack 
of 11 GRB hosts is < 7 X 10^ Mo for a 70 Myr old pop- 
ulation, and < 3 X 10^ for the maximal age. We also 
consider the previous 5 long GRB hosts at z ^ 3 (Table [2]), 
including two detections with maximal masses of 1.4 x 10^^ 
Mo rGRBQ6Q 51QB: Chary eta l. 2007) and 6.7 x 10^^ Mo 
(GRB 080607 l iChen et al.ll2010h in our M*-Z analysis. 
In comparison to these values, the typical stellar mass of the 
I GRB host sample is about 1.2 x 10^ Mo, similar to our 
stack limit. Similarly, the most massive GRB hosts at z ~ 1 
have masses of ~ 10^^ Mo, similar to our detected hosts. We 
reach a similar conclusion in comparison to the LBG sam- 
ple at z ^3: the typ ical stell ar masses of the S pitzer detected 
LBGs (iReddv et al . 2006; M agdis et al.ll201Ql) are about 10^^ 
Mo, although some of these galaxies (~ 5%) have stellar 
masses in excess of 10^^ Mo- This is similar to the distri- 
bution of our detected hosts and the 5 literature hosts. Our 
stack (maximal) limit falls below the typical stellar mass de- 
tection limit of the LBG sample. However, in the absence of a 
detailed mass function, it is difficult to estimate how the flux 
limit associated with the LBG selection compares to our de- 
tected fraction. Regardless of the exact answer, it is clear that 
deeper observations of the GRB host sample will probe lower 
mass systems than available with the LBG sample. 

4. THE MASS-METALLICITY RELATION AT Z 3 - 5 

We now turn to the primary investigation of this paper — 
the M*-Z relation at z ~ 3 - 5. In Figure \T0\ we present the 
absorption line metallicities plotted versus the stellar masses 
inferred from our Spitzer observations ("our sample"). Also 
included are the 5 GRB hosts from previous targeted obser- 
vations ("literature sample", TableS]). Of the 18 GRBs in our 
sample, six have determined [S/H] values, five have no metal- 
licity information, and six have upper limits on their metallic- 
ity from non-detections of S II as well as lower limits based 
on Si II or Si IV detections. For GRB 050908, the metalHc- 
ity upper limit (Z < 10^ Zo) is not meaningful and we only 
report a lower limit based on a C II detection. For the lit- 
erature sample, the spectrum of GRB 080607 exhibits a sat- 
urated S II line, leading to a lower limit on the metallicity; 
GRB 060223A has both upper and lower limits on the metal- 
licity based on a S II non-detection and a Si II detection, re- 
spectively; GRBs 0605 lOB and 050904 have measured metal- 
licities; and GRB 060522 has no metallicity information. 

Using these values we find a wide range of metallicities^ 
(Z ^ 0.01 - 1 .5 Zo) for the 5/7te^r-detected GRB hosts, which 
have stellar masses of ~ 2 x 10^^ Mo- This range indi- 
cates that at least some of the hosts have metallicities that 
are typical of z ~ 1-2 galaxies in the same mass range. We 
note that this range is significantly larger than the scatter in 
metallicity obser ved at low redshift, w hich is about 0.4 dex 
at logM* - 10 (iTremonti et al.l l2004 . Since GRBs probe 
the metallicities of their host galaxies along random lines of 
sight (whereas direct galaxy spectroscopic observations yield 
luminosity- weighted metallicities), this larger scatter may be 
indicative of the intrinsic scatter in the metallicities of indi- 
vidual star-forming regions in z ^ 3 galaxies. We return to 
this point in ^ 

To search for an M^-Z relation, we divide the GRBs with 
available metallicity information (either a metallicity detec- 
tion or a bounded range) into two mass bins — the 3.6 /im 

^ To transform the sulfur a nd oxygen abundances we use the Solar values 
Hsted in lAsplund et al.l (I2005h . 



detections with M* - 2 x 10^^ Mo (Group 1: GRBs 060926, 
050319, 060707, 060210, 0605 lOB) and the objects included 
in the stack (Group 2: GRBs 060526, 030323, 061 HOB, 
060115, 060906, 050730, 060206). For each group, we per- 
form a Monte Carlo simulation to estimate the mean metal- 
licity. We represent metallicity detections by Gaussian ran- 
dom variables with a mean equal to the detected metallicity 
and variance equal to the observed uncertainty, and objects 
with metallicity ranges by uniform distributions. The simula- 
tions yield nearly-Gaussian distributions for the mean metal- 
licity of objects in both bins, with (Zi) = -1.01 ±0.17 and 
(Z2)=-1.52±0.12, where the quoted uncertainties are 1 a er- 
rors on the mean. The mean metallicities of the two groups 
are measurably different at the level of about 1.8 a. 

For Group 1, the mean of the maximum inferred stellar 
masses is 4.3 x 10^^ Mo, while that of the masses inferred 
from the 70 Myr populations is 9.8 x 10^ Mo- To obtain mass 
estimates for Group 2, we scale our stack limit obtained for 1 1 
non-detections by ^/ll/l. Using the mean maximum mass- 
to-light ratio of the objects in Group 2 (3.7 x 10^^ Mo//iJy) 
yields an upper limit on the mean stellar mass of these seven 
objects of 3.7 X 10^ Mq, while the mean mass-to-light ratio at 
70 Myr (9.4 x 10^ Mo//iJy) yields a mass limit of 9.4 x 10^ 
Mq . The resulting mass ranges together with the correspond- 
ing la and 2a metallicity ranges for both groups are indicated 
by hatched regions in Figure (TO] We find that the mean metal- 
licity decreases as a function of stellar mass, an initial indica- 
tion of an M^-Z relation. We note that our averaging of the 
individual metallicities at a fixed stellar mass is similar to the 
approach taken by Erb et al. (2006) for their z ^ 23 sample 
for which they constructed composite spectra in various mass 
bins (i.e., they averaged the spectra, while we average the in- 
dividual metallicities). 

To compare our measurements with the observed M^-Z 
relations at lower redshifts, we need to ensure the use of 
a common stellar IMF and calibration of the spectral in- 
dices u sed to measur e the metallicities. We use the re- 
sults of Maiolino etal." (2008) who re-calibrated the z 0.07 
relation of Kewlev & Ellison (l2008h . the z ~ 1 relation of 
Savaglio et al. (2005), and the z ^ 2.3 relation of lErbet all 
(2006) to the Salpeter I MF. In FigurefTOl we pl ot the resulting 
M^-Z relations given in lMaiolino et al.l (|20Q8|) . which are of 
the form: 

Z = [0/H] = -0.0864 (logM* - logMo)' + /To - (0/H)o , (3) 

where Mq and Kq are the parameters of the log-parabolic fit 
to the re-calibrated data, and (0/H)o = 8.66 is the solar oxy- 
gen abundance (lAsplund et al.ll2005|) . We also include in Fig- 
ure [TOl theJVl^Zrek at z ~ 3.1-3.5 
(Maiol ino et al.ll2008l: iMannucci et al.ll2009l). alon g with the 
mean M^-Z points for z ^ 3.1 ('Mannucc i et al ][2009) and 
z ~ 3.5 (Maiolino et al. 2008). We find that our two points fall 
below the observed relations at z ^ 3.5, providing evidence 
that the galaxy M*-Z relation continues to evolve at z ~ 3 - 5, 
with our stack range probing a somewhat lower mass scale 
than the LBG studies at z 3. 1 - 3.5. 

5. DISCUSSION AND CONCLUSIONS 

We present the first study of the galaxy mass-metallicity re- 
lation at redshifts of z ~ 3 - 5 using GRB afterglow absorption 
metallicities and Spitzer follow-up observations. Five of the 
20 GRB hosts in our sample are detected above a 3 a flux den- 
sity threshold of 0.25 /iJy, corresponding to a typical stellar 
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mass of ~ 2 X 10^^ Mq. We further place a limit of < 3 x 10^ 
M0 on the non-detected hosts based on a stacking analysis. 

The rest-frame optical luminosities and derived masses are 
generally similar to those found for GRB hosts at lower red- 
shifts, but are larger than for LAEs at similar redshifts. The 
comparison to the LBG population is less certain. No GRB 
hosts in our sample are detected with Mab ^ -23 mag, while 
about 20% of the LBG sample are more luminous than this 
value. On the other hand, integration of the z ~ 3 optical 
luminosity function suggests that we expect only ~ 2 GRB 
hosts brighter than this limit in our sample, statistically con- 
sistent with zero detections. The recent detection of the host 
galaxy of GRB 080607 with a luminosity and mass similar 
to the most massive LBGs supports t he possibility t hat GRB 
hosts and LBGs are not dissimilar (C hen et al .1120 10 ). 

We find a wide dispersion in the metallicities of the host 
galaxies (inferred mainly from S II) at a fixed stellar mass 
of ~ 2 X 10^^ M0. The mean metallicity at this mass scale 
is about 0.1 Z©. The mean metallicity associated with 7 
of the 1 1 non-detected hosts, which have an upper limit of 
< 3.7 X 10^ M0, is Z < 0.03 Z©. Thus, there appears to be an 
overall decline in metallicity with decreasing stellar mass, a 
hint of an M^-Z relation. Furthermore, our two points on the 
M^-Z relation lie below the relations at lower redshifts, sug- 
gesting that the relation continues to evolve at least to z ~ 4. 
Clearly, additional observations are required to confirm and 
increase the statistical significance of this result. A sample of 
20 additional GRBs at z ^ 3 from 2007 through the present is 
available for study. This will allow us to double the existing 
sample. 

The observed range in metallicities at ~ 2 x 10^^ 
appears to be larger than the observed scatter in metallicities 
at similar stellar masses in the nearby universe. While it is 
possible that this is a real effect, we caution that this may be 
an observational artifact; GRBs probe individual sight lines 
through their host galaxies, whereas traditional methods inte- 
grate the spectrum over an aperture or slit, thereby averaging 
over many individual H II regions (weighted by their luminos- 
ity). Since we divide galaxies into two groups and compute 
their mean metallicities, in effect we achieve a similar result 
as integrating over the many H II regions in individual galax- 
ies. 

We end the paper with some cautionary notes and future 
prospects. One possible source of systematic uncertainty 
in metallicities as probed by GRBs lies in the radial abun- 
dance profile of the host galaxies. The Milky Way and M33 
display a strong abu ndance gradient (~ -0.07 dex kpc~^: 
iRolleston etaP (20001: ICionil.2009) . Similar abundance gra- 
dients (~ -0.05 dex kpc~^) h ave been found for H II regions 
in nearby spir al galaxies fe.g. JVila-Co stas & Edmunds" 1992': 
|vanZeeetaDri998). At higher redshift, Jones etal. (2010) 
find a gradient of -0.3 dex kpc~^ in a lensed system at z = 2.0; 
while this gradient is large in absolute terms, they clarify that 
it is similar to gradients in nearby spirals when the evolution 
of the effective radius out to z ~ 2 in taken into account. On 
the other hand, the LMC and SMC, which may be more rep- 
resentative of GRB hosts, display almost no radial metallic- 
ity gradient (Cioni 2009) . Since GRBs probe an unknown 
line of sight through their hosts, a strong metallicity gradient 
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combined with a preferred location for the progenitors may 
lead to a systematic bias in the resulting metallicity measure- 
ments. However, a distribution of GRBs that uniformly sam- 
ples H II regions within their hosts, coupled with the potential 
that low mass galaxies at high redshift have weak gradients, 
will negate such a bias. 

We note that a similar effect may exist in direct galaxy 
metallicity measurements. This is simply because the mea- 
sured metallicity is effectively a luminosity-weighted value, 
which therefore depends on the combined radial distribution 
and luminosities of H II regions. For example, if H II regions 
in the outskirts of LBGs were more luminous, an abundance 
gradient would lead to a biased metallicity value. Thus, as- 
suming that galaxies of similar masses have similar abun- 
dance profiles, we would expect metallicities determined by 
GRBs as an ensemble for a given galaxy mass to be represen- 
tative of the typical galaxy metallicity at that mass. We con- 
clude that the effect of metallicity gradients is minimal when 
comparing samples as a function of galaxy mass. 

A second concern is the relative calibration of absorption 
metallicities (using mainly the sulfur abundance) and nebu- 
lar emission line metallicities (using mainly the oxygen abun- 
dance). At present, the uncertainty in the solar abundance of 
these two elements (primarily oxygen) leads to at least ~ 0.1 
dex uncertainty in the relative calibration. Beyond this prob- 
lem, an additional concern is that while to first order GRB 
absorption spectra and nebular lines both trace regions of star 
formation, it is unclear how the luminosity-weighted nebular 
metallicities relate in detail to the line of sight GRB metal- 
licities (even in the absence of metallicity gradients). Thus, 
cross-calibration of the metallicities using direct spectroscopy 
of GRB hosts is of the utmost importance. This is missing at 
the present. 

Looking beyond additional Spitzer observations of existing 
GRB hosts, which from our work appear to have a detection 
yield of ~ 25% at z ^ 3, the James Webb Space Telescope 
(JWST) will provide a much deeper view of the M^-Z rela- 
tion at high redshift. For instance, the NIRCam instrument on 
JWST employing the 3.6 jim wide filter will be able to detect 
point sources at a 5 a flux of about lOnJy in a similar integra- 
tion time to our existing observations. This will allow us to 
detect GRB hosts at z ~ 3 down to a mass of ~ 10^ M©, and 
at z 6 to 3 X 10^ M0. At these limits, we should be able 
to detect the bulk of the hosts individually if they are simi- 
lar to GRB hosts at z ~ 1 . Equally important, the NIRSpec 
instrument (1-5 jim) will allow us to determine emission- 
line metallicities for some of the hosts, and hence to cross- 
calibrate the afterglow absorption metallicities. Thus, with 
on-going afterglow absorption metallicity measurements, the 
GRB sample will continue to play a key role in our study of 
high redshift galaxies. 
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TABLE 1 

Spitzer OBSERVATIONS AND INFERRED PROPERTIES OF GRB HOST GALAXIES AT Z ~ 3-5 IN THIS STUDY 



ORB Z F^(3.6/im) Pec ^ MAB(780nm) (M/F^)70Myr (M/F^Uax MvOMyr Mr 









(mag) 


(lO^^Mo/ziJv) 


(lO^^Mo/ziJv) 


(lO^^Mo) 


(lO^^Mn) 




060607 


3.075 


<0.24 


— 


> -20.14 


0.69 


2.97 


<0.17 


<0.71 


020124 


3.198 


<0.26 


— 


> -20.29 


0.73 


3.09 


<0.19 


<0.80 


060926 


3.206 


1.65 ±0.07 


4.3% 


-22.30 ±0.07 


0.73 


3.10 


1.21 ±0.05 


5.12±0.22 


060526 


3.221 


<0.25 


— 


> -20.25 


0.74 


3.12 


<0.19 


<0.78 


050319 


3.240 


0.80 ±0.09 


7.7% 


-21.53 ±0.20 


0.75 


3.14 


0.60 ±0.07 


2.51 ±0.29 


050908 


3.344 


<0.25 




> -20.31 


0.78 


3.27 


<0.20 


<0.82 


030323 


3.372 


<0.21 




> -20.14 


0.79 


3.32 


<0.17 


<0.70 


060707 


3.425 


1.10±0.10 


6.8% 


-21.96±0.16 


0.81 


3.40 


0.89 ±0.08 


3.74 ±0.34 


061 HOB 


3.433 


<0.25 




> -20.35 


0.82 


3.41 


<0.20 


<0.85 


060115 


3.533 


<0.31 




> -20.63 


0.85 


3.56 


<0.26 


< 1.10 


060906 


3.686 


<0.28 




> -20.58 


0.90 


3.76 


<0.25 


< 1.05 


060605 


3.773 


<0.24 




> -20.45 


0.92 


3.88 


<0.22 


<0.93 


050502 


3.793 


<0.21 




> -20.32 


0.93 


3.90 


<0.20 


<0.82 


060210 


3.913 


1.41 ±0.10 


5.9% 


-22.42 ±0.11 


0.97 


4.08 


1.37±0.10 


5.75 ±0.41 


050730 


3.968 


<0.26 




> -20.62 


0.99 


4.15 


<0.26 


< 1.08 


060206 


4.048 


<0.21 




> -20.41 


1.02 


4.29 


<0.21 


<0.90 


060927 


5.464 


<0.21 




> -20.83 


1.44 


6.98 


<0.30 


< 1.47 


050814 


5.772 


0.55 ±0.06 


7.6% 


-21.94 ±0.20 


1.54 


7.75 


0.85 ±0.09 


4.26 ±0.47 



Probability of chance coincidence. 
2 This is a photometric redshift dCurran et al.ll2008l) . 
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TABLE 2 

Spitzer OBSERVATIONS AND INFERRED PROPERTIES OF GRB HOST GALAXIES FROM PREVIOUS STUDIES 



ORB 


z 




MAB(780nm) 


(M/F^)70Myr 




^VOMyr 








(MJy) 


(mag) 


(lOiOM0//iJy) 


(lOiOM0//iJy) 


(IOIOMo) 


(lOi^Mo) 


080607 ^ 


3.036 


22.9±2.1 


-25.07 ±0.09 


0.67 


2.94 


15.4±1.4 


67.3 ±6.2 


060223 A b 


4.406 


<0.30 


> -20.91 


1.13 


4.89 


<0.34 


< 1.47 


060510Bb 


4.942 


0.23 ±0.04 


-20.78 ±0.17 


1.29 


5.86 


0.30 ±0.05 


1.35 ±0.23 


060522b 


5.110 


<0.20 


> -20.68 


1.34 


6.18 


<0.27 


< 1.24 


050904^ 


6.295 


<0.27 


> -21.29 


1.69 


9.12 


<0.46 


<2.46 



^ IChen eta l.'(>2010) 
b Chary et al. ( 2007) 
SBergeretaL (.20071) 



Mass-Metallicity Relation at z ~ 3 - 5 



TABLES 

Metallicities from Afterglow Absorption Spectroscopy for GRB Host Galaxies at 







z ~ 3- 


-5 IN THIS Study 








z 


lOg(A'Hl) 


on 


IVCbL VVdVclcllgul 




LZ./11J 






(cm~^) 




(A) 


(cm~^) 




060607 


3.075 


16.95 ±0.03'' 










020124 


3.198 


91 70 + 20^ 










060926 


3.206 




S II 


1253.8 


< 15.97^ 


< -1.77 








Si II 


1526.7 


> 14.56^ 


> -3.55 


UOUJZO 


^ 991 


ZU.UU ZC u. 1 J 


Q TT 


2 

vanous 


14 -1- n 9^ c 

l'4-.Z)o ZC U.ZJ 


n ^7 -U n 9S c 
u. J / zc u.zj 


Q 

\jJ\jDVy 




90 on -1- n 90 a 


TT 


1 9S'^ s 
izj j.o 


^ 1 A 1 < a 
^ 10. 1 J 


^ n 1 1 








Si IV 


1402.8 


\ 1 /I A7 a 
> 14.0/ 


\ 1 7/1 

> —1 . /4 


\JJ\jy\jO 


3.344 


i / .\J\J ZC U. iU 


CII 


1334.5 


> 13.85 ^ 


^ —0.14 


\JD\JDAj) 


D.D 1 L 


9 1 on -U n 07 d 
LL.yyj zt u.u / 


S II 


1253.8 


1 J.o^- zt U. 


1 90 -1- n 90 

—1 .zu zt u.zu 


yjyjyj lyj 1 




91 on -1- n 90 ^ 

zi .uu ziz u.zu 


S II 


1250.6 


1 A an -1- n 9 ^ 

iO. J)U ZIZ u.z 


n 1 -1- n 98 

U. iO ZIZ U.ZO 


061 1 lOR 


3.433 


22 35 + 10^ 


S II 


1253.8 


< 15.65 ^ 


< -1.84 








Sill 


1304.4 


> 15.04^ 


> -2.82 


060115 


3.533 


21.50±0.10^ 


SII 


1253.8 


< 16.12^ 


< -0.52 








Sill 


1526.7 


> 14.72^ 


> -2.29 


060906 


3.686 


21.85±0.10^ 


SII 


1253.8 


< 15.63^ 


<-1.36 








Sill 


1526.7 


> 14.25^ 


>-3.11 


060605 


3.773 


18.90 ±0.40^ 










050502 


3.793 












060210 


3.913 


21.55±0.10^ 


SII 


1253.8 


15.80±0.10f 


-0.89 ±0.14 


050730 


3.968 


22.10±0.10^ 


SII 


1253.8 


15.11 ±0.04g 


-2.13±0.11 


060206 


4.048 


20.85 ±0.10^ 


SII 


various^ 


15.21 ± 0.03 h 


-0.78 ±0.1 


060927 


5.464 


22.50±0.15^ 


SII 


1253.8 


< 16.90^ 


<-1.55 








Sill 


1260.4 


> 13.99^ 


> -4.02 


050814 


5.773 













^ Solar abundances are from Asplu nd et al.l (120051) . The metallicities have been re-derived from the quoted metal 
ion column densities, as necessary. 

^ Derived by simultaneous least-squares fitting of the S II A1250.6, A1253.8, and A1259.5 transitions. 



^ This is a photometric redshift (iCurran et al.ll2008h . 
^ Fvnbo et al. ( 2009) 
^ Hiorth et al. ( 200J) 

Thone et al. ( 20® 

^ Vreeswiik et al. (120041) 
^ Ferrero et al. (2009|) 
^ This work 
§ Ledoux et al. ( 2(23) 
^ . Fvnbo etal. (20()^ 
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TABLE 4 

Metallicities from Afterglow Absorption Spectroscopy for GRB Host 
Galaxies from Previous Studies 



GRB 


z 


log(A^Hi) 
(cm-2) 


Ion 


Rest Wavelength 
A 


log(Mon) 
(cm-2) 


[Z/H]i 


080607 


3.036 


22.70 ±0.04^ 


SII 


1250.6 


> 16.3 


>-1.5 


060223A 


4.406 


21.60±0.10b 


SII 


1253.8 < 15.3b 


<-1.4 










Si II 


1304.4 


- 15.3b 


>-1.8 


0605 lOB 


4.942 


21.30±0.10b 


SII 


1250.6, 1253.8 


15.6±0.lb 


-0.85 ±0.15 


060522 


5.110 


21.00±0.30b 










050904 


6.295 


?^21.6^ 


SII 


1253.8 


15.6±0.15^ 


_1 14+0.14 



the quoted metal ion column densities, as necessary. 

^ Prochaska et al. ( 2009j 

b Chary et al. (2007) 

^Totani et al. ( 2006) 

^ IKawaietal1f2^ 
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TABLE 5 

Field Galaxy Comparison Samples 



Sample 






z 




Detections 


Limits 


^AB> ^ 


(Mab> ' 


Depth 
















(mag) 


(mag) 


(mag) 


Reddv et al. (2006) 


2.29 


< 


z < 3.00 


2.94 


32 


6 


23.2 


-22.2 


24.9 


^eddv et al. (2006) 


3.00 


< 


z < 3.66 


3.22 


31 


5 


23.7 


-21.9 


24.9 


^agdis et al. (2010) 


2.34 


< 


z < 3.00 


2.93 


19 


3 


23.4 


-22.0 


25.5 


Magdis et al. (2010) 


3.00 


< 


z < 3.45 


3.20 


30 


6 


23.6 


-22.0 


25.5 


ShaDlev et al. (2005) 


1.48 


< 


z < 2.90 


2.30 


72 





22.8 


-22.1 


R < 25.5 


Ono et al. (2010) 


3 


•1, 


3.7^ 


3.40 


11 


261 


23.7 


-22.1 


24.8 


This work 


3.0 


< 


z<5.8 


3.43 


5 


13 


23.8 


-22.0 


24.8 b 



^ Me dian apparent magnitudes in the Spitzer 3.6 jim band dReddv et al.ll2006l : lMagdis et "ZIIMol i lOnoern] 
[20TQ ) and /T-band ( Shaplev et al. 2005i) . 
^ Corresponding median absolute magnitudes. 

^ This work on LAEs uses two narrow-band filters tuned to Lya at z ~ 3.1 and z ~ 3.7. 
^ Median apparent magnitude of our 3.6 non-detections. 
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Fig. 2. — Spitzer images of the GRB host galaxies with 3.6 12m detections. The left-hand panels show the processed images, while the right-hand panels include 
subtractions of nearby sources using GALFIT (when performed). The circles (1 radius) mark the afterglow positions. All images have the same orientation (North 
is up and East is to the left) and scale (16" on a side) with 0^^4 square pixels. 
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Fig . 5. — The luminosity distribution fo r GRB hosts at z ~ 3 - 5 from our study (black) and the 5 GRB hosts at z > 3 from previous studies (bluei lBerger et al.l 
12007: "Chary et al.ll2007l : 'Chen et al."201Q), compared with GRB hosts at z ~ 1 (red). All samples have been K-corrected to 780 nm in the rest frame. Non- 
detections (3cr upper limits) are shown as open histograms. Also shown is the stack limit at a median redshift of z ~ 3.4. The host of GRB 080607 at Mygo nm = 
-25.1 was targeted due to the large extinction inferred from the optical afterglow (Chen et al.ii201Q) . 
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Fig. 6.— Norm alize d luminosity distribu tion for the GRB hosts in our sample (black: solid=detections; dashed=limits) compared with LBGs at z ^ 3. 
IReddv et al.l f200^ and'Magdis et al." (^20T^ use Spitzer 3.6 /im observations of GOODS-N and probe deeper than our study (Table |5), while Shaplev et al.J 
( 120051) rely on ground-based K-bmd photometry. K-corrections for the SED shape (last term in equation |2) have not been applied (although the relative 
difference in K-correction should be minor and will not modify the shape of the distribution). Non-detections have been removed from the comparison samples 
for clarity. The upper panel shows the corresponding distributions of the rest frame wavelengths for each sample. 
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-18 -19 -20 -21 -22 -23 

M (AB mag) 



Fig. 7. — Same as Figure|6] but compared with Spitzer 3.6 /im observations of Lyman-alpha emitters at z ~ 3.1 and ~ 3.7 (red; lOno et al.ll201Cih . Only 4% 
of the LAE sample was detected individually (red histogram), while stacks at z ~ 3.1 and ~ 3.7 revealed much lower typical luminosities of Mab ~ -18.3 mag 
and ^ -19.4 mag, respectively. K-corrections for the SED shape (last term in equation |2) have not been applied (although the relative difference in K-correction 
should be minor and will not modify the shape of the distribution). 
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Age(Gyr) 



Fig. 8. — Upper panel: mass-to-light ratio in the rest-frame /-band from the Maraston (2005) simple stellar population models in Solar units. The models 
assume an instantaneous burst of star formation (r = 0). Lower panel: ratio of stellar mass to observed 3.6 /im flux density for (r = 0) SSP models as a function 
of age at different redshifts from z = 3 (bottom curve) to z = 6 (top curve) in steps of 5z = 0.5. The models used in this paper are for a Salpeter IMF with a red 
horizontal branch morphology and a metallicity of 0.022© . Each curve is truncated at a value that corresponds to the age of the universe at that redshift. 
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Fig. 9. — Stellar masses plotted as a function of redshift for our sample (black). Bars indicate detections with the upper and lower ends defined by the maximal 
mass-to-light ratio and a 70 Myr old stellar population, respectively. Individual limits (triangles) are plotted for the maximal masses, while the dashed curve 
indicates the 3a limit (0.25 /iJy) for a 70 Myr old stellar population. The stack limit for 11 non-detections at z ~ 3.4 is designated by the green bar (70 Myr 
old population). Also shown are the 5 previous Spitzer observations at z > 3 from the literature (blue symbols), and low redshift GRB hosts from the study 
of Leibler & Berger (2010). The additional c omparison samples include star forming galaxies at z ~ 2 (orange; Shapley et al. 2005) and LBGs at z ~ 3 (red; 
iReddv et al.ii2006 : >Maiolino et al.i.2 008: Mann ucci et alJI200^ : iMagdis et al.ii2010) . The Spitzer-detQCted GRB hosts at z ~ 3-5 have similar masses to the most 
massive GRB hosts at z ~ 1 and to the LBGs at z ~ 3. The stack limit is similar to the typical masses of GRB hosts at z ~ 1. 
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Fig. 10. — Stellar mass plotted as a function of ISM metallicity for our sample and the 5 previously observed hosts (black=detections; gray=limits). Confirmed 
metallicities (and one metallicity lower limit) are indicated by filled symbols, while metallicity ranges are shown by dashed vertical lines with open symbols. 
The red hatched regions designate la and 2a intervals for estimates of the mean metallicity at two mass bins of ~ 2 x 10^^ Mq (3.6 /im detections) and 
< 3.7 X 10^ M0 (scaled stack limit - see text). The red dashed vertical line indicates the upper limit on the mean mass of the stack for a 70 Myr p opulation. 
These data are consistent with a decline in metallicity with lower stellar mass — an M^-Z rel ation. Also shown are the relations for z ~ 0.07 ( Kewlev & EUisonl 
[2008), z ~ 0.7 ( Savaglio et al. 2005), z ^ 



^2.3 (Erb et al. 2006), and z ~ 3.1-3.5 jMaiolino et al. 2008; Mannucci et al. 2009, filled squares); the relation at z < 2.3 

are the re-calibrated values bv lMaiolino et al.l ( l2008>) . Our two data regions at z ~ 3 - 5 fall below these relations suggesting that the M*-Z relation continues to 
evolve to z ~ 4. 



